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The Clausius relation between entropy change and heat flux has previously been used to derive 
Einstein's field equations as an equation of state. In that derivation the entropy is proportional to 
the area of a local causal horizon, and the heat is the energy flux across the horizon, defined relative 
to an approximate boost Killing vector. We examine here whether a similar derivation can be given 
for extensions beyond Einstein gravity to include higher derivative and higher curvature terms. We 
review previous proposals which, in our opinion, are problematic or incomplete. Refining one of 
these, we assume that the horizon entropy depends on an approximate local Killing vector in a way 
that mimics the diffeomorphism Noether charge that yields the entropy of a stationary black hole. 
We show how this can be made to work if various restrictions are imposed on the nature of the 
horizon slices and the approximate Killing vector. Also, an integrability condition on the assumed 
horizon entropy density must hold. This can yield field equations of a Lagrangian constructed 
algebraically from the metric and Riemann tensor, but appears unlikely to allow for derivatives of 
curvature in the Lagrangian. 
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I. INTRODUCTION 

The notion of horizon entropy and thermodynamics 
was first discovered for black holes in general relativity 
(GR) and quickly generalized to other sorts of horizons. 
The origin of this thermodynamic behavior can be traced 
to local physics, and in a sense arises from the nature of 
the vacuum. This led to the observation that the Ein- 
stein equation can be derived as an equation of state for 
local causal horizons in the neighborhood of a point 'p' 
in spacetime, by imposing the Clausius relation between 
their entropy change and the energy flux across them pQ . 
In this paper we examine approaches to generalizing this 
equation of state derivation to allow for higher derivative 
contributions to the entropy and field equations. The 
Einstein-Hilbert Lagrangian is only the lowest order term 
(other than the cosmological constant) in a derivative ex- 
pansion of generally covariant actions for a metric theory, 
and the presence of higher derivative terms is presumably 
inevitable. Several approaches to including higher deriva- 
tive terms in the entropy and equation of state have been 
tried. In our view none have fully succeeded, except in 
the case of C(R) theories. Those theories are special how- 
ever, since they are trivially related to general relativity 
coupled to a scalar field by a field dependent conformal 
rescaling of the metric. Here we will explain the prob- 
lems that arise with previous proposals, and propose a 
solution that adopts aspects of some of the proposals. 

The solution differs in several ways from the original 
derivation for GR, among which are: (i) the entropy is 
compared on two horizon slices that share a common 
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boundary, (ii) the bifurcation surface lies to the past of 
the terminal point p and, (iii) the entropy depends on the 
approximate Killing vector. In particular, it has the same 
dependence on the approximate Killing vector as have the 
Noether charges associated with a Lagrangian, in analogy 
with the Wald entropy [2| |3j for stationary black holes. 
Such dependence will be referred to as "Noetheresque" , 
but by itself does not make the entropy a Noether charge. 
It makes some sense that the entropy depends on the ap- 
proximate Killing vector, because the latter determines 
the notion of stationarity and enters the definition of the 
heat flux. However, at present we can offer no statistical 
interpretation for this form of the entropy. 

The equation of state we can derive is consistent with 
local energy-momentum conservation only if the leading 
order term in the entropy satisfies an integrability condi- 
tion. This condition is satisfied if the entropy arises from 
variation of a generally covariant function with respect 
to curvature. In other words, the entropy coincides with 
a Noether charge associated with a (particular type of) 
gravitational Lagrangian. The need for such an integra- 
bility condition was anticipated in Ref. [T], since it was 
known that the entropy of stationary black hole horizons 
has this form[3H5]. We have not been able to ascertain 
whether this is the only way to satisfy the integrability 
condition, and in particular whether field equations for 
Lagrangians involving derivatives of curvature can be ob- 
tained. 

Besides the lack of a statistical interpretation, and the 
dependence of the entropy on the local Killing vector, 
a strange feature of this approach is that in the case of 
GR the entropy on a general horizon slice differs from 
the area at the same order as the relevant area changes 
in the Clausius relation. So this approach seems to lose 
contact with some of the original statistical motivation 
for the local Clausius relation. On the other hand, it is 
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quite analogous to the first law of black hole mechanics in 
generalized gravity theories. Therefore it is not clear to 
us whether this approach is purely formal, or maintains 
some significance as true thermodynamics of spacetime. 

II. THE EINSTEIN EQUATION OF STATE 

We begin with a review of the derivation of the Einstein 
equation as an equation of state, which emerges from 
black hole thermodynamics as follows. General relativ- 
ity and quantum field theory in a black hole background 
imply the so-called "first law of black hole thermodynam- 
ics" [BHE], 

dM -ft H dJ = T H dS B H, (1) 

where M is the black hole mass, the angular velocity 
of the horizon, J the angular momentum, Th = 11k/2tt 
with surface gravity n is the Hawking temperature, and 
Sbh = A/(4hG) with horizon area A is the Bekenstein- 
Hawking entropy. This relation can be viewed as a com- 
parison between two stationary black holes, but it also 
holds for small, slowly time-dependent changes of a sin- 
gle black hole. Its validity in that setting hinges on the 
fact that the evolution of horizon area is governed by 
spacetime curvature, which in turn is linked via Einstein's 
equation to the energy flux. Specifically, it relies on the 
relation R a bk a k h = 8irGT a bk a k b , where R ab is the Ricci 
tensor, T a b is the energy momentum tensor of matter, 
and k a is a 4-vector tangent to the horizon generating 
null geodesies. 

The term "first law" is actually a misnomer for ([I]) . In 
thermodynamics that name refers to energy conservation, 
dU = SQ + SW, where dll is the internal energy change, 
SQ is the heat flow into and SW is the work done on the 
system. Instead, the thermodynamic nature of is the 
Clausius relation 

SQ = T5S (2) 

between the heat flow and the entropy change. In the 
black hole context, energy that flows across the horizon 
is, in effect, "heat" , since after crossing the horizon its mi- 
croscopic nature is effaced for an outside observer. Since 
the horizon is a causal barrier, it is a "perfect dissipa- 
tor" [9J. The heat flux is 

5Q = dM-Q H dJ= f {-T abX a )dH\ (3) 

where the integral is over the horizon and x a = df+^lnd^ 
is the horizon-generating Killing vector. 1 That is, the 



We (unfortunately) use spacetime signature ( — h***+). The 
horizon integration measure is dHi, = —k^dVdA, where V is 
the afHne parameter along the horizon generators, k b = (dy) 
is tangent to the generators, and dA is the area element of a 
constant V horizon slice. 



heat is the energy flux conjugate to the spacetime sym- 
metry that translates along the horizon generators. The 
work term SW in the usual first law of thermodynamics 
has no analog in the relation ([!]). Although "first law" is 
not an appropriate name, we will use it here since it is 
entirely standard terminology. 

Close to a black hole horizon the Hawking tempera- 
ture becomes the Unruh temperature ha/2ir for station- 
ary, uniformly accelerated observers, and the Clausius 
relation takes on a local form (how local depends on how 
slowly the changes take place) whose validity requires 
that the Einstein equation hold, as mentioned above. 
Moreover, there is good reason to believe this applies not 
just to black holes but to any causal horizon. (See e.g. 
|10) for a review of the arguments, and for further 
discussion and clarification.) 

The idea of Ref. [T] (see also [T^] for a slight reformu- 
lation) was that, conversely, the Einstein equation can 
be derived, as an equation of state, by requiring that 
the Clausius relation hold for the entropy of sufficiently 
small patches of all local causal horizons (LCHs) in space- 
time. Any such horizon H is defined as the boundary of 
the past of a patch of (D — 2) -dimensional spacelike sur- 
face, where D is the spacetime dimension. (LCHs will 
be defined precisely in what follows; for now it suffices to 
remain somewhat vague.) In that derivation the heat is 
taken as the boost energy flux, defined with respect to 
an approximate local boost Killing vector field £° as 

SQ= f {-T ab e)dH\ (4) 

and the temperature is taken as the Unruh boost tem- 
perature h/2ir. This is natural since the Minkowski vac- 
uum state of quantum fields is thermal with respect to 
the boost Hamiltonian at this temperature, and any state 
looks like the Minkowski vacuum at sufficiently short dis- 
tances. The entropy is taken to be A/Iq, where Iq is some 
UV length scale. (A compelling case can be made that 
the origin of horizon entropy is quantum entanglement 
of vacuum correlations across the horizon. Since this is 
dominated by UV degrees of freedom, yet finite, dimen- 
sional analysis suggests that the leading order contribu- 
tion should scale in this way.) The change in horizon area 
in a small neighborhood of a point p can be related to the 
Ricci tensor via the Raychaudhuri equation. If the hori- 
zon is chosen so its expansion and shear vanish at p, and 
£ a is chosen so it too vanishes at p, the Clausius relation 
is then seen to require T a bk a k b — (h/(2Trl^))R a bk a k b at 
every point in spacetime and for all null vectors k a . To- 
gether with energy conservation V a T a b = this implies 
the Einstein equation 

Rab — \Rg a b — A.g a b = 8irGoT a b, (5) 

where the value of Newton's constant is determined by 
the entropy density 1/Zq to be Go = Iq/4H, and A is an 
undetermined cosmological constant. 
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III. PREVIOUS APPROACHES TO 
INCLUDING HIGHER DERIVATIVES 

There have been a number of attempts to include 
higher derivative terms in the equation of state de- 
rived from causal horizon entropy; all these attempts 
can be divided into two broad classes: The first deals 
with a specific theory, L(R) gravity (see Refs. [T2HT1] ). 
while the second class studies a more general case(e.g. 
Refs. [13 [16]). The L(R) gravity, where the Lagrangian 
depends only on the Ricci scalar, is equivalent to gen- 
eral relativity with an auxiliary scalar field JTTJ [TB] and 
is therefore the simplest generalization possible. Even 
in that case, the derivation of the field equation from 
thermodynamics of local horizon meets considerable ob- 
stacles. 

Ref. [H] considered the case when the horizon entropy 
density is proportional to an arbitrary function f(R) of 
the spacetime Ricci scalar R. As was done in the deriva- 
tion of the Einstein equation, the approximate Killing 
vector £ is chosen to vanish at a terminal point p on 
the horizon. The Clausius relation equates the entropy 
change to the heat flux SQ, which is the flux of boost 
energy current — T ah £& across the horizon. Since £ van- 
ishes at the terminal point, the rate of change of entropy 
with respect to the affine parameter must also vanish 
there. However, at a generic spacetime point p the gra- 
dient V Q i? is non-vanishing, so the horizon entropy will 
have a nonvanishing rate of change unless the change of 
R is balanced by a change of horizon area. Thus it is 
necessary to adjust the terminal surface of the horizon 
so that the expansion 9 of the horizon generators at the 
equilibrium point has the nonzero value 9 p = — ///. The 
Clausius relation involves a 9 p term via the Raychaudhuri 
equation, and one ends up with a field equation that is 
inconsistent with matter energy conservation unless this 
term is identified as internal entropy production due to 
a bulk viscosity proportional to f(R) and added to the 
entropy balance law. Hence, instead of equilibrium ther- 
modynamics, a non-equilibrium approach is needed. The 
resulting field equation is the one that follows from the 
Lagrangian L, where / = dL/dR. That is, L is the La- 
grangian for which f(R) is the Wald entropy [3]. (The 
properties of Wald entropy will be briefly reviewed in 
section 
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In Ref. |13j , this entropy production term is interpreted 
as a separate contribution to the heat flux from the ad- 
ditional scalar degree of freedom present in L{R) gravity. 
Then it is possible to derive the field equations using 
only reversible thermodynamics as in the case of general 
relativity. Although such an interpretation may offer a 
possible understanding of entropy production terms for 
the specific case of L(R) gravity, it is unclear how to 
generalize that for a broader class of theories. 

In Ref. [14] it was proposed that the need for internal 
entropy production could be eliminated by adopting 
the instantaneous boost invariant (IBI) prescription 
for dynamical horizon entropy proposed in [5J. This 



prescription goes as follows. In the neighborhood of a 
spatial slice £ of a causal horizon one constructs, by 
a unique recipe of dropping selected terms in a Taylor 
expansion of the true metric in an adapted coordinate 
system, a new spacetime metric that has an exact 
boost Killing field for which E is a fixed point set. The 
entropy density of the slice £ was taken in 14J to be 
a scalar formed from the IBI metric associated to S. 
In particular the case with entropy density f(R) was 
discussed, but it is not clear to us that all contributions 
to the change in this entropy were taken into account. 
The difficulty arises because the IBI metric changes 
with the horizon slice, so the R on each slice is defined 
with respect to a different metric. Also, even if the 
field equation arrived at in [14] were correct, it refers to 
the curvature of the IBI metric, and it remains unclear 
how this is related to the curvature of the original metric. 

Next we would like to discuss attempts to obtain the 
field equation from local horizon thermodynamics beyond 
L(R) gravity. Ref. [TS] starts from the Wald entropy for- 
mula 3] for a stationary black hole, but applied to a 
local horizon. A formula for the entropy of a [(D — 1)- 
dimensional] patch of a local horizon is written as a sur- 
face integral over the boundary of the patch. It is unclear 
to us why an entropy should be assigned to a (D — 1)- 
patch rather than only to a (D — 2)-slice of the horizon. 
Moreover, in the case of a stationary black hole horizon 
this integral vanishes, whereas the entropy of the horizon 
surely does not vanish. Hence this quantity cannot be in- 
terpreted as the entropy of the horizon. The variation of 
this entropy is then considered, and a formula is written 
in terms of an integral of a covariant directional deriva- 
tive of part of the integrand of the previously mentioned 
integral. It is unclear to us why this expression should be 
interpreted as the change of entropy. A similar formula 
is written for the boost energy flux across the horizon. 
Also, the binormal on horizon cross sections is identified 
with the covariant derivative of the approximate Killing 
vector. Whereas this can certainly be satisfied to lowest 
order, given the small changes involved in the Clausius re- 
lation it would be necessary to check whether neglecting 
the differences at subleading order is justified. Another 
issue is that the approximate Killing vector £ a in these 
calculations is treated as if it satisfies the Killing identity 
V a V(,£ c = R d abc £,d exactly. Although this is not possible 
in general, we will show in the present paper that the 
identity can be satisfied to the required order, provided 
one restricts to a narrow neighborhood of a particular 
horizon generator. 

In Ref. [THj the proposal is to define the entropy of 
horizon slices as the integral of a Noether potential 
formally identical to a potential introduced in [3] , except 
that in the former case no matter fields are included in 
the Lagrangian with which the potential is associated. 
Using Stokes' theorem, the entropy change between two 
slices of the horizon is then expressed as an integral 
of the correspondig Noether current over the enclosed 
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horizon patch. However, since slices of a LCH have 
boundaries, the entropy change also involves a surface 
integral over the null outer boundary, unless the two 
slices have their boundary in common. This contribution 
was missed in Ref. |16j . In addition, as in Ref. |15j . 
the Killing identity was used without justification at 
an order it cannot be expected to hold in a generic 
spacetime. Thus this derivation too is not complete. 
The present paper starts from a similar but more general 
form for the entropy, and fills these gaps in the argu- 
ment by studying slices that do have their boundary in 
common and thus enclose a compact patch of the horizon. 

Finally we note that Refs. [HI |2Qj (see also [21] and 
references therein) discuss a thermodynamic interpreta- 
tion of gravitational field equations beyond general rel- 
ativity in terms of a local entropy balance law. A mat- 
ter entropy flux across the horizon is associated with a 
small spatial volume and is related to the boost energy by 
the Clausius relation. This is then equated to a gravita- 
tional entropy change in the volume, which is constructed 
from the Noether current associated with a Lagrangian 
L[g a b, R a bcd\- This balance law at a point is then shown 
to imply that the gravitational field equations hold, thus 
giving a thermodynamical interpretation of those equa- 
tions. The approach we will pursue in this paper uses 
similar ingredients, but we have different objectives. We 
start by assigning a horizon entropy functional to a small 
but finite sized horizon slice, and investigate what prop- 
erties it must have if it is to satisfy the Clausius relation 
with the matter boost energy flux. 



IV. IMPOSSIBILITY OF GENERALIZING THE 
NON-EQUILIBRIUM APPROACH 

A natural question to ask is whether the non- 
equilibrium approach of Ref. |12| can be extended be- 
yond theories for which the horizon entropy density s is 
a function only of the Ricci scalar. In this section it will 
be shown that the non-equilibrium approach of Ref. [12] 
cannot be extended beyond the results found there. For 
this discussion, we adopt the geometric set up of [TJ 
and we assume that the entropy associated with the LCH 
takes the form 



s dA, 



(6) 



where s is some arbitrary scalar function and the integra- 
tion is over a (D — 2)-dimcnsional spacelike slice of the 
horizon with 'area' element dA. The change in entropy 
from one horizon slice to another is 



SS = 



dX 



6s dX dA, 



(7) 



where A is an affine parameter on the local horizon gen- 
erators and 9 = d(ln dA) /dA is the expansion of the gen- 
erators. 



The Clausius relation asserts that the entropy change 
SS is equal to SQ/T, where T = h/2n is the boost tem- 
perature, SQ is the heat flux Q through the horizon and 
£ a is the approximate Killing vector which vanishes at 
the final equilibrium point p. Since the heat integrand 
vanishes at p, the Clausius relation can hold only if the 
SS integrand also vanishes there. This requires that the 
expansion at p is nonzero, to wit 



1 ds 
s dA 



(8) 



In the case of GR s is a constant, so this condition states 
that the horizon must have vanishing expansion at p. The 
Raychaudhuri equation can be used to find 0(A). If the 
horizon shear does not vanish at p we must add internal 
entropy production due to shear viscosity to the Clausius 
relation |12j . Assuming the shear at p vanishes, we have 
6 = e p + A(-0 2 p /(D - 2) - R ab k a k b ) + 0(A 2 ), where fc° is 
the affine horizon tangent vector (d/dA) a , and A p = 0. 

Now we require that the Clausius relation hold for all 
local horizons through p, i.e. for all k a . If the entropy 
density is independent of k a , as would be the case if it is 
a spacetime scalar, then we obtain the tensorial equation 

D — 1 i 2tt 

sR ab - V Q V b S + — -S S, a S, b +*.9afc = — Tab, (9) 

L> — 2 rj 

where \& is some scalar function. This function can be de- 
termined by imposing energy conservation, i.e. V a T Q t> = 
0, which yields the condition 



~^sR, a +d a Ds - I — \ s 1 s, a s, b 



(10) 



Since ^ a is the gradient of a scalar, a solution exists only 
if the right hand side is also the gradient of a scalar. 

In the case of GR s is a constant, so \& = —sR/2 is the 
unique solution up to a constant (the cosmological con- 
stant). If s = s(R) is a function only of the Ricci scalar, 
then the first term on the right hand side of ( 10 ) is a gra- 
dient, but the last term is not. In the non-equilibrium 
approach of Ref. [T2] the Clausius relation is replaced by 
an entropy balance relation SS = SQ/T + SSi, where SSi 
is an internal entropy production term due to bulk viscos- 
ity, which is proportional to the square of the expansion 
and which cancels the last term in (10 1. This yields an 



equation of state that coincides with the field equations 
for L(R) gravity. For more general spacetime scalar en- 
tropy densities, e.g. s(R a b c d), even the first term is not 
a gradient, so there is apparently no way to satisfy an 
entropy balance law, even allowing for internal entropy 
production. Hence, the methods devolved in [T] and [12] 
can not be generalized beyond L(R). 

The situation is even more problematic if the entropy 
density is not a spacetime scalar but, for example, is 
constructed from the intrinsic curvature of the (D — 2)- 
dimensional horizon slice, as in Lovelock gravities [22j . 
In that case, we do not even have a tensorial relation- 
ship like Eq.([9|. The Raychaudhuri equation then does 
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by this choice. The role in the new approach will be 



explained in Sec. VIII 



Figure 1: Causal spacetime diagram of the local causal hori- 
zon (LCH) and Killing field in the old (left) and new (right) 
setups. Each point in the diagram represents a patch of space- 
like 2-surface. The boundary of the past of a patch of 2- 
surface through p defines the LCH. The arrows indicate the 
flow lines of the local Killing field £ a , which is of third order 
on the bifurcation surface B, and which vanishes at p (left) or 
po (right). The central horizon generator T runs from po to 
p. The "heat" is the boost energy flux across the horizon. In 
the new setup the Killing vector is timelike within the horizon 
and spacelike beyond, while the old setup is the opposite. 



not seem to be of any help at all, so a very different ap- 
proach is needed. In this paper we will present such an 
alternative approach, which involves an entropy density 
that depends explicitly on the Killing vector, a different 
choice of horizon slices, and the application of Stokes' 
theorem. 



The role of the approximate Killing field £ a in the ap- 
proach of Ref. \1[ is to define the heat in the Clausius 
relation Q. The bifurcation surface of this Killing field 
was previously taken to coincide with the future bound- 
ary through p, whereas in the first law of global horizon 
mechanics the bifurcation surface lies to the past of the 
perturbation. Moreover, the approximate Killing field 
was spacelike rather than timelike in the region outside 
the horizon, so it was necessary to think of the heat as 
going into the reservoir behind the horizon (see Fig. [lj. 3 
However, this reservoir is not observable on the outside, 
so should play no role in the outside thermodynamics. 
It seems much more satisfactory to place the bifurcation 
surface to the past, so that the Killing field is timclikc 
outside the horizon, and the reservoir can be thought of 
in direct analogy with the thermal atmosphere of a black 
hole, below the stretched horizon. 4 

It appears at first that modifying the location of the 
bifurcation point of the Killing vector will change the 
heat flux and, in the original approach of Ref. [T], ruin 
the derivation of the Einstein equation from the Clausius 
relation. However, that is not what happens. The old 
Killing vector on the generator through p (see Fig. [IJ 
was ££ ld = — Afc;,, where A is the affine parameter that 
vanishes at p. The new Killing vector vanishes instead at 
Po and is given by = (A — \o)kb, where Ao is the value 
of A at po- The corresponding boost energy currents are 
thus 



NEW CHOICE OF APPROXIMATE 
KILLING VECTOR 



T ab£old = X T ab k b 



-T a % = (X Q ~\)T ab k b 



(11) 
(12) 



It turns out that the new approach to the thermody- 
namic derivation of field equations, even in the case of 
general relativity, calls for a small but essential adjust- 
ment in the choice of approximate Killing vector. This 
adjustment amounts to locating the bifurcation surface 2 
at a time earlier than the terminal point, rather than be- 
ing coincident with the terminal point as originally for- 
mulated in Ref. [1 . The change also makes the applica- 
tion of the Clausius relation more closely analogous to 
the first law of black hole mechanics [35] and resolves 
an uncomfortable aspect of the earlier derivation. In 
this section we explain the thermodynamic motivation 
for the new choice of approximate Killing vector, and 
show how the original derivation of Ref. [T] is modified 



Although these differ, their integrals from the bifurcation 
point to p are the same, since 



AdA 



Xl/2 



(A - A) dX. 



(13) 



With the heat defined using £°, applying the Clausius 
relation to the horizon interval [Ao, 0] in the limit Ao — > 
thus yields the field equations as described above. This 
corresponds to a transition between a stationary state at 
A = Ao (where £ a vanishes) and one at A = (where the 
expansion and shear vanish). In analogy to the physical 
process version of the first law of black hole mechanics, 
the Clausius relation holds only when applied to this en- 
tire interval. 



2 In the current context, the term "bifurcation surface" is used 
somewhat loosely, mainly to fix attention to the analogous role 
played by the true bifurcation surface of a quasi-stationary black 
hole in the physical process version of the first law |23| . What 
we mean is a slice on which the components of the approximate 
Kiling vector are third order in an appropriate coordinate sys- 
tem (see Sec. | VI P\ . Accordingly, the approximate Killing vector 
vanishes at the center of the bifurcation surface (see Fig. 111. 



3 We use the term "outside" to refer to the region accessible to 
observers to the past of p, like the observers outside a black hole 
horizon. For a cosmological horizon, the standard terminology 
would be opposite: outside a cosmological horizon is the region 
that can not be seen. 

4 This same point was made recently in Ref. |24| . with reference 
to the choice of the "observer" who is defining the heat flux. 
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Note that in the approach of Ref. [lj, with the new 
choice of Killing vector, one must still impose the condi- 
tion that the expansion 6 vanishes at p, because that con- 
dition is used in applying the Raychaudhuri equation to 
obtain the area change. By contrast, in the Noetheresque 
approach of the present paper the expansion at p will turn 
out to play no role whatsoever. 

VI. PROPERTIES OF THE LOCAL HORIZON 
AND KILLING VECTOR 

In this section we spell out the detailed geometric 
construction and properties of the local causal horizon 
(LCH) and Killing vector. 

A. Local causal horizon 

We define a local causal horizon H as follows. Consider 
any spacetimc point p in a D-dimensional spacetime, and 
let E p be any small patch of spacelike (D — 2)-surface 
through p. The boundary of the past of E p in the neigh- 
borhood of p has two components, each of which is a null 
surface generated by a congruence of null geodesies or- 
thogonal to E p . The LCH H is defined as one of these 
components. 

As mentioned in the previous section, it will not mat- 
ter for the present approach whether the congruence is 
stationary at p. However, if p is to be a stationary point, 
the expansion and shear of this congruence must vanish 
there. That is, the extrinsic curvature of E p must vanish 
at p, which is equivalent to saying that E p is generated 
by geodesies at p. For concreteness we will go further and 
assume that E p is fully generated by geodesies emanating 
from p. This will be convenient for studying concrete ex- 
amples of the approximate Killing vector and of entropy 
values. 



B. Null normal coordinates 

To establish the existence of an approximate Killing 
vector with the required properties, we will employ a 
"null normal coordinate" (NNC) system [35] adapted to 
the LCH. This is an explicit realisation of the coordi- 
nate systems first introduced in and is defined 
as follows. At a spacetime point p an orthonormal set of 
(.D— 2) spacelike vectors {e^}, A = 1,2, .., (D-2), is cho- 
sen, and a (D — 2)-surface E p is generated by geodesies 
with tangent vectors at p in the space spanned by {e^}. 
The point reached on such a geodesic at unit affine pa- 
rameter is assigned the coordinates x A , when x A e a A is the 
tangent vector in that paramctrisation at p. This defines 
standard Riemann normal coordinates on E p based at p. 
A pair of future pointing null vector fields (k a , l a ) normal 
to E p is chosen on E p , normalized such that k a l a = — 1. 
Each point r in a small enough spacetime neighborhood 



of p lies on a unique geodesic orthogonal to E p at some 
point q. Let the tangent to that geodesic at q be given 
by Vk a + Ul a when r lies at unit affine parameter from 
q. The NNCs of r are then defined by 5 

x? = (U,V,xf). (14) 

These coordinates are defined uniquely up to a rotation 
of the (D — 2)-frame at p and a g-dependent rescaling 
of k a , with inverse rescaling of l a . The horizon H is 
the surface U — 0, restricted to V < 0. In particular, 
the central horizon generator T is the coordinate curve 
U = 0, x A = 0, V < 0, with bifurcation point p e T 
at V = Vq < 0. Note that on H the coordinate V is an 
affine parameter along the null generators. The details of 
the construction, including the derivation of metric coeffi- 
cients and Christoffel symbols, are presented in Ref. [25] . 
It is also shown there that the ambiguity in the choice of 
k a can be exploited to make the coordinates locally in- 
ertial at p and to further specialize the properties of the 
metric components. A preferred choice is made by start- 
ing with the Z?-dimensional Riemann normal coordinates 
at p and adjusting them. This also induces a particular 
choice of affine parameter V. 

C. Horizon slices 

The Clausius relation applied to a horizon refers to the 
entropy change SS between two times. Those times are 
spacelike hypersurfaces, one to the future of the other, 
which intersect the horizon in two slices. For a black 
hole with a compact horizon, the two slices may bound 
a cylindrical region of the horizon. 

For a local causal horizon, the considered process must 
be local, since the LCH is not even well defined except in 
a small neighborhood of the terminal point p. To local- 
ize the process we can restrict attention to cases where, 
rather than pushing the spacelike hypersurface forward 
in time everywhere, it is deformed to the future only in 
a small neighborhood of the bifurcation point pq. 6 Then 
the two corresponding horizon slices also coincide every- 
where except in a small region. Parts of two such slices, 
Eo and E, are depicted in Fig. [2} E corresponds to the 
V = Vq surface and E lies to the future. If we truncate 
the horizon slices outside the region where they differ, 
their union E U E forms the closed boundary of a patch 
of the horizon. This will allow the difference of the en- 
tropies on the two slices to be computed using Stokes' 
theorem. 



5 We use greek indices, and upper case latin indices, for coordi- 
nates and coordinate components in the NNC coordinate sys- 
tem. Lower case latin indices are reserved for abstract indices on 
tensors. 

6 Hypersurface deformations of this sort were considered in 
Ref. |27| in the context of proving a local form of the Gener- 
alized Second Law. 
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v = o 




Po(r-o) 



Figure 2: The local causal horizon is part of the boundary of 
the past of E p , a spacelike D — 2 dimensional surface. The 
vertical lines are horizon generators with affine parameter V 
and tangent vector k a , and F is the central generator. E p is 
geodesic at p, so the expansion 9 and shear vanish at p. The 
horizontal lines represent constant V slices. The local Killing 
vector vanishes at the bifurcation point po- The entropy is 
compared on two slices, Eo and E, which differ only in a 
compact region. 



In fact, we shall need to further restrict the choice of 
horizon cuts so that they bound a narrow region of the 
horizon, because only then can the Killing identity be sat- 
isfied to sufficient accuracy for the approximations made 
in the derivation of the equation of state to be valid 7 . By 
"narrow" , we mean that the ratio of the width in x A to 
the length in V goes to zero in the limit as po approaches 
P- 



D. Local Killing vector 

Next we define precisely the approximate Killing vector 
£ a that plays a central role in defining both the heat flux 
and, in the Noetheresque approach of Sec. |VIIIB[ the 
entropy density. We will refer to this vector field as the 
"local Killing vector" . 

Of course a general curved spacetime has no Killing 
vectors. Nevertheless, in a small enough neighborhood 
of any point, any spacetime is approximately flat. In 
particular, in local inertial coordinates at p, the metric 
components take the form g a p — rj a p + 0(x 2 /L 2 ), where 
i] a p is the Minkowski metric, x denotes the coordinates, 
and L characterizes the shortest radius of curvature of 
the spacetime at p. (Here and below we use greek indices 



7 An alternative option would be to restrict to a horizon patch 
that is symmetric in the transverse coordinates, to cancel con- 
tributions from terms linear in those coordinates. However, we 
consider this symmetry restriction to be artificial. 



to refer to components in a particular coordinate sys- 
tem, reserving lower case latin indices for abstract ten- 
sor indices.) An approximate boost generator £ a can be 
defined in terms of local inertial coordinates by the for- 
mula for an exact fiat spacetime boost generator, e.g. 
(x, t, 0, 0) in Minkowski coordinates. This vector satisfies 
the Killing equation 



(15) 



exactly at p, and to 0(x) near p, but in general V( Q £b) 
will have 0{x 2 ) terms. This is adequate for our purposes. 

However, our application of the Clausius relation is 
also sensitive to the 0(x) part of the Killing identity 



V a Vf,£ c — R abc^d- 



(16) 



For a true Killing field, this identity follows from the 
Killing equation (15 1. Conversely, the Killing identity 



implies the Killing equation if the latter holds at one 
point, because of the antisymmetry of the Riemann ten- 
sor in the last index pair. Our computations will rely on 
this identity being satisfied at O(x), but for the approxi- 
mate Killing field defined above it will generally not hold 
at that order. We can try to modify the definition of this 
vector field so as to satisfy both the Killing equation and 
the Killing identity at 0(x), but in fact this is not possi- 
ble. However, narrowing our sights, £ b can be chosen so 
that the Killing identity holds to this order (in fact it can 
be chosen to hold exactly) on the single horizon generator 
r that ends on the terminal point p, and this turns out 
to be good enough. In effect, the local Killing symmetry 
can be extended away from p to a better approximation 
along a single null generator V than across the whole 
LCH. This calls to mind null Fermi coordinates [28], but 
we have used the NNCs to describe the situation since 
they are better adapted to the LCH. 

Given that we confine the horizon to a narrow region 
surrounding the central horizon generator T, the integrals 
appearing in the Clausius relation will be dominated by 
their integrands evaluated on T. As such, any conditions 
that £ a may need to satisfy in order for the Clausius 
relation to lead to a consistent equation of state will be 
conditions imposed on T. 

As motivated in Sec. [Vj the local Killing vector £ a is 
taken to have a bifurcation surface Eo — or at least a bi- 
furcation point — at po to the past of p, where it vanishes 
and where its covariant derivative generates boosts in the 
plane orthogonal to the bifurcation surface. The NNCs 
of po are (U, V, x A ) = (0, Vb,0). It will be convenient to 
shift the afflnc parameter such that its origin coincides 
with the bifurcation point, i.e. we define V = V — Vb, so 
Po lies at V = 0. We also want £ a to be approximately 
tangent to the generator T that connects po to p, so that 
the evolution of the LCH can be interpreted as a small 
perturbation of a stationary background, which justifies 



the use of the Clausius relation. 8 More specifically, the 
equation of state derivation will require that, at least to 
lowest order, £ a coincides with Vk a , as would be the case 
on a Killing horizon. And, finally, the derivation will re- 
quire the Killing equation and the Killing identity to hold 
at 0(V). Expressed in NNC components, the full set of 
requirements that must be imposed on £ a is: 



{k a lp - l a kp)\ 



V( Q £/3)| r 

V Q V^ 7 | r = (R s aP ^ s )\ r + 0(V 2 ). 



V5 V + 0(V 2 ), 
0(V 2 ), 



(17) 
(18) 
(19) 
(20) 



As alluded to above, £ a may be subjected to rather 
stronger conditions under which it approximates a true 
Killing vector more closely. In particular, we may choose 
£ a such that the Killing identity (16) holds exactly on 



r. This can be demonstrated by a perturbative argu- 
ment which is given at the end of this section (see also 
[25] for further details). The identity then implies that 
the Killing equation also is exact on T, provided it holds 
exactly at one point. Another consequence is that £ a 
becomes exactly tangential to T, as will be shown next. 
Hence our clai m i s t hat all terms of quadratic or higher 
order in V in (18) - (20) may actually be set to zero. 



If a vector field ( a satisfies the Killing identity along a 
geodesic with affine tangent vector v a , it follows that 



V a V a KVfeCc) = RdabcCv a V b 



(21) 



Eqn.(21) is an ordinary differential equation of second 



order and has a unique solution along the geodesic, once 
£ a and v b VbC a are given at an initial point. It is easily 
verified that sv a is a solution to (21), where s is any 



affine parameter on the geodesic, and v Vb(sv a ) 



sv 



and u b V h C a 



at one point, it 



Hence if ( a = 
follows that £ a = sv a everywhere on the geodesic. Now 
recall that k a is an affine tangent to T, and the local 
Killing vector is chosen to vanish at po and to satisfy 
(kJb ~ lah)L„- Taking v a = k a thus yields 



Ml 



V a 6 

v b Vb£, a \ po — v a \ po , so we may conclude that £ a = Vk a 

everywhere on T. 9 

As regards the Killing equation, we recall that the 
NNC system is local inertial at p. As discussed at the 
beginning of the current section, it is then possible for 
the Killing equation (15) to hold at 0(x), and this turns 



out to be consistent with the Killing identity on T. In 
NNCs, the approximate Killing equation implies that the 
components of the local Killing vector are of the form 
= VS V - U5%+0(x 3 ). 

In summary, the main properties exhibited by our 



choice of local Killing vector are given by (17) together 
with 



Tlr 



V5y, 

0(x 2 ), 



(22) 
(23) 
(24) 



To conclude this section, we turn to the perturbative 
argument that the Killing identity ( 16 1 can be satisfied 
exactly on T. We specify £ a in a neighborhood of V by 
its (covariant) components in NNCs as the Taylor series 

£ Q = USl - VS% + C^ a x^SP + D^ Sa x p x^x s + . . (25) 

where x a = x a — VoSy and in particular V — V—Vo. Sim- 
ilarly, on r the deviation away from the Killing identity 
is written as a power series in the affine parameter V. As 
shown in [35], the latter series may be set to zero order 
by order through an appropriate choice of the expansion 
coefficients occurring in ([25]) . 



More specifically, at the linear order (20) needed for 



the equation of state derivation, one finds 

Capf ~ Vq x (Riemann components at p), (26) 
Dvafi-y ~ (Riemann components at p). (27) 

Thus the terms that are quadratic and cubic in x 
in the expansion ( 25 ) both contribute to £ Q only at 
0(x 3 ), counting Vq as 0(x). This is consistent with the 



Killing equation (23) at O(x), from which one finds that 



quadratic order terms must be absent in ( 25 ) 



For the purpose of the general derivation of field 
equations all we need to know is that the conditions 
(171-pO]) can be met. For the purpose of computing the 



actual entropy (which will have a dependence on the 
local Killing vector, see Sec. VIII B ) and comparing with 



the area in the GR case or with other expressions in 
higher derivative gravity, the detailed form of the local 
Killing vector is generally needed. For more explicit 
expressions of the higher order coefficients in ( 25 ) we 
refer to E>]. 



8 Moreover, if we were to let the local Killing vector and the tan- 
gent to the generator vary independently, it appears the Clausius 
relation would be overly restrictive. 

9 It will generally be inconsistent with the Killing identity on T to 
choose § a to be tangential to all the generators of a small horizon 
patch. Moreover, if the horizon expansion does not vanish at 
the bifurcation point, such a choice would also be inconsistent 
with the Killing equation imposed at linear order in a spacetime 
neighborhood. 



VII. BLACK HOLE ENTROPY AS NOETHER 
CHARGE 

In this section we review Wald's expression for black 
hole entropy in terms of the Noether charge[H |3]. This 
expression will motivate the form of LCH entropy that 
we adopt. 

Consider a diffeomorphism-invariant Lagrangian field 
theory of gravity in arbitrary dimensions. A variation of 



9 



the action is always equal to a sum of terms proportional 
to the field equations and a surface term. When the 
variation is induced by a vector field £ a , this equality 
can be expressed as the statement that a certain current 
j a is closed, i.e. V a j a = 0. The current is closed even 
when evaluated off-shell and takes the form 



■a =Q a _ LS a _ fa gab _ T ab\ ^ + 



(28) 



where the symplectic current 9 a stems from the sur- 
face term, L is the Lagrangian scalar in the action, 
2E ab - T ab = is the metric field equation and the 
dots indicate off-shell terms corresponding to the mat- 
ter field equations. The Noether current J a is defined as 
ja _ Qa _ L^a^ anc j on g^g]] coincides with the closed 

current j a . Furthermore, since the on-shell Noether cur- 
rent is closed for any possible vector field £ a , it can 
be expressed as J a — 2V&Q a6 , where the antisymmet- 
ric tensor Q ab is a local function of the fields and their 
derivatives [29], The tensor Q ab is referred to as a Noether 
potential. 10 

Since the vector field £ a only enters the expression for 
a Noether potential when taking Lie derivatives of tensor 
fields with respect to £ a , any Noether potential is of the 
form 



Q ab [£] = W abc £ c + P abcd V c £ d 



(29) 



apart from the possible addition of a total divergence. 
More specifically, 

Q ab [£] = W abc £, c + X abcd V [c £ d] + Y ab + W c Z abc , (30) 

where the tensors W, X, Y, Z are locally constructed from 
the dynamical fields, Y is linear in Lie derivatives with 
respect to £, and Z abc is totally antisymmetric. The de- 
composition of Q in terms of W, X, Y, Z is not unique, 
and in addition the Noether potential has three sources 
of ambiguity, coming from the freedom to add a total 
divergence to the Lagrangian, the symplectic potential, 
or the Noether potential itself. Using all this freedom Y 
and Z can be set to zero, X can be chosen [3] as 



X 



abed 



dL 



dRabcd 

and W is then given by [30] 



(31) 



W abc = 2W d X abcd + matter terms + (32) 

where the dots indicate terms that stem from derivatives 
of the Riemann tensor in the Lagrangian. 



10 In this paper we use tensor notation. The formalism in Refs. [2ll3l 
makes use of differential forms, which are dual to these tensors. 
The relation between a p-form F and the corresponding tensor 
F is F ai .,.a p = F bl "' bn ~ p <^b 1 ...b n _ p a l ...a p , where e is the volume 
element. The tensor corresponding to the exterior derivative dF 
is (n - p)V a F b ^- b «-P- ia . 



If the theory admits stationary black hole solutions 
with a regular bifurcation surface, variations away from 
these solutions satisfy the first law (1), with the entropy 
Sbh defined by [3J 



2tt 



1BH 



Q ab [x]NabdA. 



(33) 



Here, £ denotes any slice of the stationary horizon, N a b = 



2k 



[ a lh] is its binormal (normalized as N a f,N 



ah 



-2), 



dA is the (D — 2) dimensional area element, and x a is 
the horizon-generating Killing vector normalized to unit 
surface gravity. The integral is referred as the Noether 
charge and is invariant under the three sources of am- 
biguity [31 [S] . Invoking the stationary symmetry, it was 
further shown in these references that the entropy may 
equally be expressed by substituting Q ab [x] in (33 1 with 



X abcd N cd . Thus the W term does not contribute to the 
black hole entropy. For general relativity, Sbh becomes 
the familiar Bekenstein-Hawking entropy A/ AUG. 



VIII. CLAUSIUS RELATION AND EQUATION 
OF STATE 



In this section we adopt the assumption that LCH's 
have an entropy of Noetheresque form, i.e. the entropy 
density depends on the local Killing vector in the same 
way that the Noether potential ( 29 ) depends on the hori- 
zon generating Killing field. With this entropy, we shall 
find that the Clausius relation applied to all LCHs, to- 
gether with the local conservation law for the matter 
stress tensor, can be satisfied provided that i) the en- 
tropy density can be identified with (the gravitational 
part of) a Noether potential of some Lagrangian, and ii) 
the fields satisfy the metric field equation for that La- 
grangian. 11 This is the equation of state for the corre- 



sponding entropy function. 



A. Clausius relation for a local causal horizon 

This Clausius relation for a LCH is analogous to the 
first law for stationary black holes. However, the latter 
relies on the stationarity of the black hole horizon, as well 
as the fact that a horizon slice is a complete boundary 
component of a spacetime slice, neither of which hold for 
a LCH. Nevertheless, it turns out that the construction 
can be localized enough to make the Clausius relation at 
least well defined. First, as explained in Sec. |VI C| and 
Fig. [2j we compare the entropy on two LCH slices S and 
Eo that share a common boundary, so that together they 
form the boundary of a local patch H of the horizon. The 



11 Note that for the equation of state to be well defined the matter 
stress tensor need not be the functional derivative with respect 
to the metric of an action. 
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lack of stationarity looks more problematic because, for a 
dynamical black hole horizon in a general theory, there is 
no well defined notion of entropy available 12 . The prob- 
lem is that there is no Killing vector and, if a vector field 
is somehow selected, the ambiguities in the definition of 
the Noether charge will make the entropy ambiguous, un- 
like for a stationary black hole horizon. 13 On the other 
hand, locally there is always an approximate Killing vec- 
tor, the "local Killing vector" constructed in Sec. VI D 
This turns out to act enough like a Killing vector to make 
the Clausius relation well defined in a certain local limit, 
once a form for the entropy is adopted. 

In the context of horizon thermodynamics it is natural 
to expect that the entropy is an extensive quantity. Thus 
we assume that it can be expressed as an integral of a 
(D — 2)-form over a spacelike slice E of the horizon. We 
adopt the dual description, and express the entropy as 
the integral 



change integrand of (35), multiplied by T, must be equal 
to the 0(x) heat flux integrand, up to 0(x 2 ) terms. That 
is, 



(h/Tr)V b s ab k a = T Qb 6 k a + 0{x 2 ). 



(37) 



This relation is imposed at all points p and for all null 
directions at p. 

We emphasize that the limit p —¥ p is taken in a formal 
sense only, in order to identify the leading order contri- 
bution for a small region. From a physical point of view, 
it makes no sense to consider an arbitrarily small region, 
because quantum fluctuations of the metric presumably 
invalidate our semiclassical considerations at sufficiently 
short distances. 



B. Entropy density of Noetheresque form 



S 



„ab 



N ab dA, 



(34) 



where the entropy density s is an antisymmetric ten- 
sor, N ab is the binormal to the slice and dA is the area 
element on the slice. The change of this entropy between 
the two horizon slices is 



5S 



S — So 



s ab N ab dA 
sus 

2 / V b s ab dH a 

JH 

-2 f S7 b s ab k a dVdA 7 
Jh 



(35) 



where Eo is taken with the opposite orientation to E, and 
in the third line we have used Stokes' theorem (for details 
of Stokes' theorem on a null surface see Ref. |31j). 

On the right hand side of the Clausius relation, dS = 
SQ/T, we have the heat flux 



SQ= f {-T ab )£, b dH a 
Jh 



(36) 



divided by the Unruh boost temperature T = H/2tt. The 
heat flux integrand is proportional to the local Killing 
vector which vanishes at po and is thus of 0(x) in NNCs 
in the neighborhood of p®. The Clausius relation is im- 
posed in the limit po — > p, which means that the entropy 



Inspired by Refs. [TBI [TJ5], we propose a local entropy 
density of the same form as the Noether potential of Eqn. 



(291, where £ a is the local Killing vector, but a priori no 



restrictions are placed on W and P, other than that they 
are constructed locally from the dynamical fields. We em- 
phasize that, although we are using a similar notation, 
we do not assume at this stage that the divergence of 
the local entropy density is a Noether current associated 
to a Lagrangian, in contrast to the entropy density of 
stationary black holes. This form of entropy density can 
be viewed as the most general one that depends linearly 
on a local Killing vector that satisfies the Killing identity, 
and contributes to the entropy change at 0(x 2 ). As men- 
tioned in the introduction, since the local Killing vector 
determines the local notion of equilibrium and heat, it 
is not entirely unnatural that the entropy would depend 
on the Killing vector. On the other hand, this entropy is 
not strictly intrinsic to the horizon, and has no immedi- 
ate statistical interpretation that we are aware of. 
We thus assume the entropy density takes the form 



with 



ab _ ^Cfab 



(38) 



= w abc £ c + P abca V c £,d 



(39) 



In the second line of (39), the tensor P is split into an 



antisymmetric part X and a symmetric part Y: 

j£-abcd j£ab[cd] yabed ^ab(cd) 



(40) 



12 General relativity is an exception, where the area of the horizon 
slices is a good candidate for the dynamical entropy. This is 
motivated by the area theorem, as well as by the form of entan- 
glement entropy; see Sec. [TT] 

13 Although the proposed dynamical entropy of Ref. [3] (in terms 
of the instantaneous boost invariant metric) manages to bypass 
the ambiguities, it most likely does not satisfy a second law 32 . 



We could of course absorb the factor of 2w/h into the 
definition of Q ab . The only reason we factor it out here 
is so that the notation will make the analogy with the first 
law for stationary black hole horizons more transparent. 

We now show that the symmetric part does not con- 
tribute to the Clausius relation. The contribution of 
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yabcd tQ Vb gab Jg 

v b (y o6cd v (c e d) ) - V b Y abcd V (c £ d) + Y abcd W b W c U 



(41) 



Using the approximate local Killing equation ([23]), we see 
that the first term on the RHS is of 0(x 2 ). The second 
term is symmetric in cd, so without any change we may 
add a term in the Riemann tensor, yielding 



Y ahcd ( V & V c £ d 



R 



J 



bed 



if 



(42) 



which according to the approximate Killing identity ( 20 ) 
is of 0(x A ), where x A are the transverse coordinates of 
the NNC system. The contribution of (41) is therefore 



negligible in the limit of small, narrow horizon patches. 

Thus, at least for the part that contributes to the en- 
tropy change, our proposal for the entropy density re- 
duces to 2ir/h times 



gab = W at>c£ c + X abcd w ^ d 



(43) 



where W and X are, so far, unspecified tensors, con- 
structed locally from the dynamical fields, antisymmetric 
in the first two indices, and in addition X is antisymmet- 
ric in the last two indices. Note that the X term is of 
0(1), and the W term is of 0(x). Had we kept the Y 
term it would have been of 0{x 2 ). 

Next we show that the Clausius relation requires that 
W is a combination of divergences of X. The divergence 
of ( 43 ) can be written as 



V b Q Q0 = (V r W arb + X arsi R° rst )£ b 

+ X arst (V r V s tt-R b r S tZb) 

+ (W ast +V r X arst ) V.&. (44) 

Again, the first term is proportional to the local Killing 
vector and hence of 0(x), while the second term involving 
the Killing identity may be neglected on narrow patches. 
It is precisely here that the need to impose the Killing 
identity arises. Were this identity not satisfied by £ a , 
the second term would make a contribution at the same 
order as the first term, which would depend on the higher 
order coefficients in the local Killing vector expansion. 
This would ruin our extraction of the local field equation 
from the Clausius relation and would probably make it 
impossible to consistently impose the Clausius relation 
at all for such local Killing vectors. 

In the third term, the symmetric part in (st) is of 0(x 2 ) 
by virtue of the Killing equation, but a priori the an- 
tisymmetric part is of O(x ). On the other hand, the 
integrand in the heat flux is of 0(x), so in order to be 



consistent with the Clausius relation (37), the O(x ) part 
of the third term must vanish. It follows that the local 
Clausius relation can only be satisfied at all spacetime 
points if W and X are related via 



W a 



V r X° 



0. 



(45) 



Because W is antisymmetric in its first two indices, this 
relation completely determines W as a function of X, to 
wit 



W a 



= V s (X sarb + x sbra + x srba ) 



(46) 



Thus only the first term in ( 44 ) survives 



Notice that although the W term does not contribute 
to the entropy at po? n or at p in the limit po —> p, its 
rate of change is comparable to that of the X term, so 
that — in contrast to the stationary comparison version 
of the first law of black hole mechanics [3] — it makes an 
important contribution to the Clausius relation, which 
could not be satisfied without it. 14 



C. Equation of state 



Substituting ( 46 ) into the first term of ( 44 ) , and using 



the fact that £ a oc k a on the central generator, the validity 
of the Clausius relation (37) for all k a implies 



R [a rst X b ^ rst + 2\7 r \7 s X ( - a ^ b)r 



ab 



_rpab 

2 



(47) 



where $ is some scalar function that may depend on the 
metric and curvature. Note that (47) follows from the 



Clausius relation irrespective of the values of expansion 
and shear anywhere on the horizon patch. 

As was done in Sec. |IV| we impose local conservation 
of energy- momentum to determine the function Then 
Eqn. pf) leads to 



(R (a rst X b ^ st + 2V r V s X^ s ^ . (48) 



In order for such a $ to exist, the right hand side must 
be the gradient of a scalar. This integrability condition 
further constrains the nature of X, which so far is only 
required to be anti-symmetric in both the first and second 
pair of indices. 

If the integrability condition is satisfied, then the left 
hand side of (47) is a divergence free tensor constructed 
from the metric and its derivatives. One way to obtain 
such a tensor is from the variational derivative of a scalar 
action functional with respect to the metric, 5I g [g]/Sg ab , 



14 One might think that perhaps a total divergence can be added 
to the entropy density such that the W term is cancelled. That 
is, if a tensor Z abcd = Z^ ab ^ d could be found such that W abd = 
\7 c Z abcd , then we would have Q ab = (X abcd - Z abcd ) V c £d + 
V c(Z abcd ^d)- The total divergence term would not affect the 
changes in the entropy from one slice to another with common 
boundary, so for the Clausius relation this would be equivalent 
to having W = and replacing X by X — Z. The existence of 
such a Z is at first glance conceivable since, according to | |46| l, W 
must be a divergence. However, the required total antisymmetry 
of Z is not shared by the combination of X's in i |4(i[ i. Only the 
divergences of Z and X enter the equation, but still it seems 
unlikely that such a Z exists in general. 
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which is automatically divergence free. In fact, it was ar- 
gued in Ref. [33] that all such tensors arise in this way. If 
so, then the 'Clausius equation' d47| is precisely the equa- 



tion of motion that derives from the action I g 
with an undetermined cosmological constant. 



1. Relation between entropy density and Lagrangian 

We can be more specific about the relation between 
the entropy and an action whose equation of motion 



D. Examples 



L(R) theories 



We now illustrate our result by a simple example, with 
X given by 



x a. 



bed 



' ac bd 

.9 9 



ad be 

9 9 



(51) 



where / = f(R, Rij, Rijki) is an arbitrary scalar function, 
is (j47j). If a gravitational Lagrangian L[g ab , R a bcd ] is a This expression satisfies all the symmetry requirements 



scalar formed algebraically from the metric and Riemann 
tensor, the n the corresponding equation of motion is pre- 
cisely (|47|, 15 with 



j^abed 



dL 

dRabcd ' 



(49) 



and $ = L/2. Note that if X is assumed to have this 
form, it has all the symmetries of the Riemann tensor. 
Using these symmetries one can show that our condi- 
tion ( 46 1 that determines W in terms of X becomes 
exactly the same as Eqn. (32 1. Thus, with the suffi- 



cient condition (49) assumed, our conjectured entropy 



density is nothing but the specific choice discussed in 
Sec. |VII| for the Noether potential associated to Bj Lei - 
grangian L[g ab , R a bcd \. 16 

For a Noetheresque entropy density with X of the form 



(49), and W satisfying (145]) , the divergence in Eq. (44) 



is nothing but the closed current (28) of Sec. VII for the 



(50) 



Lagrangian scalar L[g ab , R a bcd ]. That is, 
2V b Q ab = a - LC - 2E ab £ b , 

where E ab is the variational derivative of L, given by the 
LHS of eqn ((4Tj) . Specifically, the first term of the RHS of 
(44) is given by —L!; a — 2E ab £ ib , while the symplectic cur- 
rent 9 a is given by the last two terms. We further note 
that the symplectic current vanishes for a Killing vec- 
tor, while it was found here to be negligible for the local 
Killing vector. Similar remarks were also made in [19]. 

Finally, we point out that this particular entropy den- 
sity does not include a dependence on matter fields. How- 
ever, such a dependence should perhaps not be ruled out 
a priori. 



15 To show this is a bit tricky. One can regard the Lagrangian 
as a function of g a i, and R a bcd j an d use the identity dL / dg a ^ = 
—2R^ a rat 8 b )dL/dR qre t- This identity can be easily established by 
exploiting the fact that the Lagrangian can also be considered as 
a function of the (2, 2) tensor R ab st , without any explicit depen- 
dence on the metric or R qrs t, and taking the partial derivatives 
using the chain rule and R t = g aq g br Rqrst- For a detailed 
derivation, see Ref. ; 34 . 

16 We note that not all Noether potentials associated to a given 
Lagrangian satisfy the condition it is curious that the local 
Clausius relation restricts the form of the Noether potential in 
this manner. Perhaps this stems from the lack of exact symmetry 
in a general spacetime background. 



imposed so far, but from Eq. (48) we obtain 



(52) 



which can not be integrated unless the function / de- 
pends only on the Ricci scalar. Thus, we find a strong 
restriction on the form of the entropy density if the Clau- 
sius relation is to be consistent with energy conservation. 
For example, the leading order term in the entropy den- 
sity can not be, say, 1 + aR ab R ab . 

If / is only a function of R, it can be written as / = 
—dL/dR for some function L(R), and then we have $ = 
L/2 (the arbitrary additive constant freedom in $ can 
be absorbed into L, and corresponds to a cosmological 
constant). With / of this form, Eq. (47) becomes 



L'(R)R ab - V a V b L'(R) 



UL\R) - ^LJ g ab = \T ab - (53) 



This is identical to the equation of motion that results 
from the Lagrangian L(R) + L ma tter- In the case of GR, 
L = R/(16ttG), and we obtain the usual Einstein equa- 
tion. 

In the case of generic L{R) theories, we have thus ob- 
tained the field equation as an equation of state from 
the Clausius relation, without the need for a term repre- 
senting internal entropy production due to bulk viscosity, 
unlike in the non-equilibrium framework of Ref. |12j . As 
described in Sec. |IV| in that approach the change of the 
entropy density f(R) coming from the gradient of the 
background R at the equilibrium point is cancelled by a 
choice of non-zero horizon expansion, which then leads 
to the bulk viscosity term in the entropy balance equa- 
tion. Instead, with the Noether charge entropy density, 
the W term is chosen so that the Clausius relation can 
be satisfied without internal entropy production. 



2. Lovelock theories 



Suppose that X abcd has no higher than second deriva- 



tives, and that also V a X abcd = 0, so (|46j) implies W abc = 
0. Then it follows that the tensor on t ic left hand side 



of the field equation (47) is second order in derivatives. 



Moreover, after imposing the integrability condition (48 1, 
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we infer that this tensor, built from the metric and its 
first and second derivatives, must be identically diver- 
gence free. The only such tensors come from the metric 
variation of a Lovelock Lagrangian|35j . so the assumed 
properties of such an X abcd , together with the Clausius 
relation, imply that X abcd arises as in (491 from a Love- 
lock Lagrangian of the form L = L[g a t,, R a h 



v = o 



bcd\ 



3. General relativity 



In example (|51| with the choice / 
46l 



implies that W — 
density obtains from Eqns 



(167TG)- 1 , Eqn. 
The corresponding entropy 



-(8fiG)- 1 Vl a ^]. The steps of Sec. 



(|43|) and (|38 | and r eads 
then 



VIIIC 



„ab 
b GR 

lead to the Einstein equation of state. 

For the Killing vector Q we have = 2fc[ Q Z b J + 

0{x 2 ). Any horizon slice has a local binormal of the form 
N a t, = 2/s[ Q /fc] + k\ a ma where m b is some spacelike vector 
tangent to the horizon and therefore orthogonal to k a . 
For generic slices, the resulting entropy (34 1 coincides to 
leading order with the area in units of 1/4HG. 

The difference in entropy between two slices requires 
a more careful analysis. On a slice of constant affine 
parameter V, the entropy will differ from the area at 
subleading order, as it does on any slice. The coefficients 
of the terms of subleading order are given in terms of 
the coefficients C and D of Eqn. (25), many of which 



are in turn set by the Killing identity ( 20 ) imposed on 
r. These coefficients can be shown to be such that the 
subleading order terms, in the case of a slice of constant 
V, involve only transverse coordinates. Therefore, when 
comparing two slices of constant affine parameter, the 
subleading terms drop out and the entropy difference to 
leading order coincides with the area difference. 

However, for the Clausius relation of Sec. |VIII A| we 
need to compare two slices with common boundary. In 
that case, on at least one of these slices the mismatch 
between entropy and area at subleading order does de- 
pend on V, in such a way that the entropy and area 
differences will not generically coincide to leading order, 
unlike in the first law of global horizon mechanics in GR 
or the original Einstein equation of state derivation pQ. 
The easiest way to see this is to employ Stokes' theorem 
([35]) with the entropy density s$ R = -(8HG)- 1 ^ a ^ . 
To evaluate the integrand of (35), we use ( f44| and recall 
that the approximate Killing vector was chosen such that 
( 44 ) is dominated by its first term. In this way we can 



express the entropy difference as 



SS « / R ab k a k b (V 

JH 



Vo) dVdA, 



(54) 



where the integration range H is a patch of the LCH 
enclosed by two slices. This integral will be equal to the 
area difference to leading order only for a very special 
choice of slices, as we now describe. 

Suppose that the first slice S is the surface V — V , 
and the shape of the second slice £ is such that the en- 



v = v„ 



Figure 3: A local causal horizon patch enclosed by a special 
pair of slices. The entropy change associated with Squ = 
— (8hG)~ 1 V' a £ 6 ' is equal at leading order to the area change 
over AhG if (i) the horizon expansion and shear vanish at p, 
and (it) the patch between the slices Eq and E nearly coincides 
with the interval [Vb,0] across the entire patch. 



closed horizon patch nearly amounts to the entire inter- 
val [Vo,0] (see Fig. [3]). Then, as indicated in (13), the 
entropy difference may cquivalcntly be approximated by 
SS s» j H R a bk a k b VdVdA. If in addition the expansion 
and shear of the horizon generators vanish at p, then the 
Raychaudhuri equation tells us that SS ~ J H 6 dVdA. 
Since the integration range is nearly the entire interval 
[Vb,0] this evaluates to SS w A(V = 0) - A(V = V ), 
where A(V) denotes the area of a slice of constant affine 
parameter V. The area of £ nearly coincides with 
A(V = 0) because its steep sides are "nearly null" and 
hence contribute negligibly to the area. Clearly, if the 
two slices are not chosen in this special fashion the en- 
tropy difference will not coincide to leading order with 
the area difference. 

Nevertheless, if the entropy difference between two 
generic slices with common boundary obeys the local 
Clausius relation, the Einstein equation must follow as a 
consequence. We can illustrate the way this works with 
a trivial example in flat spacetime. A small patch of a 
light cone may be considered to be part of a LCH, en- 
closed by two spacelike slices with common boundary. 
The boost energy flux is identically zero, so the Clau- 
sius relation tells us the entropy on both slices must be 
identical. Employing Sq R , this may be easily verified. 
However, the area of these slices is clearly not identical. 

The entropy density Sq r also shows it is not only 
thermodynamically more natural to place the bifurcation 
point of the local Killing vector field to the past of the 
terminal point (rather than coincident with the terminal 
point) as was discussed in Sec. [VJ it is actually required 
if we want the contribution to the entropy on a constant 
V slice in the GR case to correspond to the area rather 
than minus the area. 

To see why this is so, consider the fact that in the 
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local Noether charge approach we have derived the Ein- 
stein equation as an equation of state, without specify- 
ing where the bifurcation point lies. However, if we want 
the entropy to agree with the area (to leading order), 
then s ab N a i, must be a positive constant. In the case 
of GR, the entropy density is s ab oc -V[ a £ 6 ] ~ -V a £ b , 
so s ah N a i, ~ —k a l h V a £fc. Since to leading order the 
causal horizon is a Killing horizon, we have k a V a £,b oc 
±£{, oc ±fc&, where the + sign holds when the bifur- 
cation point lies to the past, so that the Killing vec- 
tor is stretching in the future direction, and the — sign 
holds when the bifurcation point lies to the future, so the 
Killing vector is shrinking (recall that we take the Killing 
vector to be future pointing on the horizon). Hence 
s ab N a f, ~ qpZ °kb — ±1. To have a positive entropy we 
must therefore take the bifurcation point in the past. We 
have no option to reverse the sign of the assumed entropy 
density, since that would lead to a held equation with the 
opposite sign for the gravitational constant. That is, a 
negative entropy change would arise from a positive heat 
flux. 

On the other hand, choosing the local Killing vector 
in this way entails a disturbing feature. If the proposed 
entropy represents the state of the horizon, one would 
expect it, and its difference, to depend only on quantities 
intrinsic to the horizon. However, the entropy difference 
is given by the boost energy flux, which depends on the 
local Killing vector, which in turn depends on a freely 
chosen location Vq of the bifurcation surface So- This is 
unsatisfying, as Vq is a freely chosen parameter and does 
not encode any intrinsic property of the LCH. 

Finally, we point out that the entropy density Sq R is 
only one among an infinite number to give rise to the Ein- 
stein equation since, by Stokes' theorem, addition to the 
entropy density of any total divergence will have no influ- 
ence on the entropy difference This is in keeping with the 
fact that the Clausius relation is a thermodynamic rela- 
tion, from which the value of the statistical entropy can 
only be deduced up to addition of an arbitrary constant. 



E. Including derivatives of curvature 

On general grounds, one may expect strong gravita- 
tional fields to be governed by a field equation derived 
from a Lagrangian that includes covariant derivatives of 
the Riemann tensor. In the present context, this raises 
the question whether such a field equation may arise as 
an equation of state corresponding to an entropy den- 
sity of the form (43 1, introduced in Sec. 



vanish for a true Killing vector, it must be expected to 
be of order one in a general spacetime. Furthermore, it 
does not appear that relation (461 can be salvaged by 



VIII B As was 



discussed in that section, in order for this to be possi- 
ble the tensors W and X must be related by Eqn. (46 1 



This may be possible, but we point out that the specific 
choice of Sec. |VII| for the Noether potential of such a La- 
grangian does not satisfy the relation ( 46 1 . The problem 



appealing to the freedom to add total divergences at the 
different stages in the construction of the Noether poten- 
tial. If this is true it would mean that the strict analogy 
to the first law of black hole mechanics breaks down, 
since that law applies for any diffeomorphism-invariant 
Lagrangian. 

However, it would still leave open the possibility that 



entropy densities of the form (43) exist that do satisfy 



relation ( 46 1 and give rise to the field equations of high 



is that the divergence of this Noether potential contains a 
Lie derivative of the Riemann tensor with respect to the 
local Killing vector. Whereas such a derivative would 



derivative Lagrangians, but that are not Noether poten- 
tials associated with such Lagrangians. As shown in the 
Appendix, the integrability condition that must be satis- 
fied in order for such entropy densities to exist (whether 
or not they are Noether potentials) may be written in 
the form of a pair of tensor equations of first order, at 
least for Lagrangians that depend on no more than first 
order covariant derivatives of the Riemann tensor. 



IX. CONCLUSION 

The question we addressed is whether gravitational 
field equations with higher derivative terms can be de- 
rived from the Clausius relation applied to a higher 
derivative horizon entropy. First we discussed prob- 
lems that arise with previous approaches to this problem. 
Then we adapted the starting point of one of those ap- 
proaches, and assumed that horizon entropy depends on 
an approximate local Killing vector in a way that mimics 
the diffeomorphism Noether charge that yields the en- 
tropy of a stationary black hole. We showed that the 
problems can all be avoided by a careful choice of the 
nature of the horizon patch to which the Clausius rela- 
tion is applied. In particular, the Clausius relation must 
refer to the change of entropy between two slices of the 
horizon that together form the complete boundary of a 
patch, and this patch must be narrow enough to neglect 
violations of the Killing identity. We exploited a power 
series expansion in a coordinate system adapted to the 
horizon to establish the required properties of the local 
Killing vector. 

Together with matter energy conservation, the Clau- 
sius relation applied to all such local horizon patches 
leads to an integrability condition on the assumed hori- 
zon entropy density. We showed that this condition can 
be satisfied if the latter is in fact a Noether potential 
associated with a Lagrangian constructed algebraically 
from the metric and Riemann tensor. In that case the 
Clausius relation implies that the field equation for that 
Lagrangian holds. We have not proved that this is the 
only way to satisfy the integrability condition, but that 
may be the case. In particular, the field equation for a 
theory with derivatives of curvature in the Lagrangian is 
unlikely to be obtained in this way using for the entropy 
density a Noether potential derived from the Lagrangian, 
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although it might conceivably arise from a different en- 
tropy density. 

The higher derivative extension of the equation of state 
derivation was achieved in this paper at the cost of intro- 
ducing a dependence of the entropy on the choice of local 
Killing vector. Whereas this dependence occurs at sub- 
leading order, the entropy difference depends on the local 
Killing vector at leading order. It is therefore not clear 
to us whether the derivation has any thermodynamic sig- 
nificance. We regard it as a positive answer to a techni- 
cal question, but its physical interpretation remains ob- 
scure. Perhaps the steps we have taken are a valid part 
of a picture that will only become clear once so far over- 
looked subtleties are taken into account. For example, 
the arbitrariness of an additive constant in thermody- 
namic entropy, the contribution from the entanglement 
entropy, and the related need to regularize by some kind 
of subtraction or comparison might play an important 
role in formulating the local thermodynamics of the vac- 
uum. Then again, it may just be that the contribution of 
higher curvature terms to a gravitational field equation 
cannot be sensibly captured at a local thermodynamic 
level. 
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W (0) 



r> -\rabcd 



a : [be] - b : [ac] - c : [ab] (A4) 



and 



a : [be] = Z a:[b def R c ^ . 



(A5) 

Due to the terms a : [be] etc., the relation (46) is not 



satisfied for Q?k . We now suppose that the entropy den- 
sity consists of the Noether potential Qfty with a further 
Noetheresque potential added. That is, we take 

s ab = 2n/hQ a \ (A6) 



where Q ab = Q$+Q$ and Q£* = W^f^+X^fV C U 
The divergence of Qff!s is given by 



(A7) 



where E ah is the variational derivative of L, the symplec- 
tic current reads 



= Z' 



:bcdt 



cde 



AV d Xf bcd 



L (o) 



A abc )V ( 6&), 



(A8) 
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Appendix A: Integrability condition for equations of 
state corresponding to Lagrangians 

L[g a t, R a bcdi V ' jR a bcd]- 



and A is a combination of terms of the form ZR. If the 



entropy density ( A6 ) is to give rise to the field equation 
derived from L, the divergence of Q?A must take the form 

OY7 /l^ rya.hcdd f r> i rpc a 

^bW(l) — L^Kbcde + f t, 



F abc V ib Z c) , 



(A9) 



by 



For notational convenience, we define the tensor Z [5] 



Z f:a 



bed 



dL 



d(y f R abcd y 



(Al) 



Then the choice of Noether potential of Sec. |VII| for a 
Lagrangian L[g ab ,R a bcd ,V f R a bcd ] is given by 



where 



\^abcd i V7 V f:abcd 

X (0) =—*E + V/Z J 

Ottabcd 



(A2) 



(A3) 



for some tensors F. From the symmetries of the Riemann 
tensor and the definition of the Lie derivative, this can 
be seen to be equivalent to the system 



W, 



a[bc 



(i) 



Vvadbc 



^ def ' 



V b W$ c 



\rabde r>c rya-.bdef \-?c r>b 

A (l) H bde- ~ Z V K def 



Fg a (AlO) 



In conclusion, if entropy densities of the desired form 
exist, there must exist an integrating function F such 
that (A10) has a solution. There is a relation between 



the integrating function F and the integrating function 
$ of Eqn. f47). Namely, comparing (|A7!-(|A9[) with (Ell) 
reveals that in this case $ = L/2 + F/2. 
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